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We investigate numerically different phases that can occur at half filling in the lowest and the 
first excited Landau levels in wide-well twodimensional electron systems exposed to a perpendic¬ 
ular magnetic field. Within a twocomponent model that takes into account only the two lowest 
electronic subbands of the quantum well, we derive a phase diagram that compares favorably with 
an experimental one by Shabani et al. [Phys. Rev. B 88, 245413 (2013)]. In addition to the 
compressible composite-fermion Fermi liquid in narrow wells with a substantial subband gap and 
the incompressible twocomponent (331) Halperin state, we identify in the lowest Landau level a 
rectangular Wigner crystal occupying the second subband. This crystal may be the origin of the 
experimentally observed insulating phase in the limit of wide wells and high electronic densities. In 
the second Landau level, the incompressible Pfafhan state, which occurs in narrow wells and large 
subband gaps, is also separated by an intermediate region from a large-well limit in which a similar 
rectangular Wigner crystal in the excited subband is the ground state, as for the lowest Landau 
level. However, the intermediate region is characterized by an incompressible state that consists of 
two four-flux Pfafhans in each of the components. 


I. INTRODUCTION 


The occurance of correlated electronic liquids and crys¬ 
talline phases in Landau levels (LLs) gives rise to a 
plethora of intriguing effects. Beyond the fractional 
quantum Hall effect,which is due to the formation of 
homogeneous incompressible liquids described with great 
success by trial wave functions,one hnds Wigner crys¬ 
tals (WCs) that are responsible for the insulating phase 
at very low hllings of the lowest LL as well as more 
exotic bubble and stripe phases,where the latter is 
responsible for an observed anisotropy in the longitudi¬ 
nal magnetoresistance. These correlated phases arise 
due to the quantization of the kinetic energy of twodi¬ 
mensional electrons in a perpendicular magnetic held B 
into highly degenerate LLs. The hlling of these LLs is 
then characterized by the hlling factor v = riei/ns, in 
terms of the electronic density riei and the density of 
states riB = eB/h in a LL. Unless an integer number 
of LLs is completely hlled (integral values of v), the 
electronic interactions remain as the only relevant en¬ 
ergy scale in the system, apart from the disorder po¬ 
tential which can at hrst be discarded in the descrip¬ 
tion of the phases. These interactions are precisely 
at the origin of the above-mentioned correlated phases, 
e.g. the fractional quantum Hall states (FQHSs) which 
occur at certain “magic” fractions of the hlling factor 
{v = 1/3,1/5, 2/5, 5/2,... and many others). In higher 
LLs, the competition between electron solids and incom¬ 
pressible liquid phases is at the origin of the reentrant 
integer quantum Hall effect,where FQHSs at different 
hlling factors are separated by insulating WC or a bubble 


crystal. 

The nature of all these ground states and the competi¬ 
tion between them are not only determined by the hlling 
factor V, but also by the precise form of the effective in¬ 
plane interaction projected onto a given LL. The form 
of the latter is affected by the index n of the partially 
hlled LL, by screening, but also by the “thickness” of the 
quantum well, which hosts the twodimensional electron 
system (2DES). This situation occurs regularly because 
the highest mobility 2DES, which allow for the resolu¬ 
tion of the more fragile FQHSs, are obtained in wide 
quantum wells with a typical width ranging from 40 to 
100 nm. For experimentally relevant magnetic helds, 
this corresponds to several (up to 20) magnetic lengths 
(.B = yJh/eB ~ 26 nm/y^i?[T]. The well width can 
even lead to the formation of novel FQHSs, such as the 
ones observed at jz = 1/2 and 1/4 in the lowest LL.^®“^^ 
The reminiscence with even-denominator states in bilayer 
quantum Hall systems^^ hints at a multi-component ori¬ 
gin later confirmed theoretically both for symmetric^^^^^ 
and asymmetric^®’^® quantum wells in terms of a (331) 
Halperin state^° that competes with the compressible 
composite-fermion Fermi liquid (CFFL)^^ and, under 
certain circumstances,with a Pfafhan.®’® Fur¬ 
thermore, these liquid phases compete with an insulating 
phase in the limit of very large well widths. 

In the present paper, we investigate the competition 
between quantum liquids and WCs at half hlling in wide 
quantum wells within a twocomponent model originally 
proposed in Ref. 23. Within this model, only the two 
lowest electronic subbands, due to the well conhnement 
potential, are taken into account whereas higher sub- 
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bands are assumed to be unoccupied. The relevant pa¬ 
rameters are the width of the quantum well and the sub¬ 
band gap that separates the lowest from the first ex¬ 
cited electronic subband. Whereas these two parame¬ 
ters are in principle not independent of each other, they 
can be tuned effectively via the magnetic held, which 
determines at the same time the magnetic length 
and the Coulomb energy e^/dTreeot's that are the natu¬ 
ral scales in the system. We investigate various relevant 
FQHSs with the help of variational Monte-Carlo calcu¬ 
lations (VMC)^^ in the spherical geometry whereas the 
energies of the electron-solid phases are calculated within 
the Hartree-Fock approximation (in the planar geome¬ 
try). In spite of the difference in these approaches and 
in the geometries, the different energies can be compared 
directly in the thermodynamic limit, which we extract 
from the usual 1 /N scaling of the VMC results. The ob¬ 
tained phase diagrams allow us to identify the different 
phases, which are responsible for the compressible, in¬ 
compressible and insulating phases identified in magneto¬ 
transport measurements at half filling in the lowest LL.^^ 
Most saliently, we obtain an insulating rectangular WC in 
the excited subband as the ground state in wide quantum 
wells with low subband gap, both in the lowest and the 
first excited LL. This phase indicates a tendency to form 
stripes at half filling and adds to the numerous different 
types of WCs observed in wide quantum wells at other 
filling factors.Whereas the intrinsic twocomponent 
state responsible for the fractional quantum Hall effect 
in the lowest LL is a (331) Halperin state, in agreement 
with previous studies,our VMC favor an unusual 
FQHS in the second (or first excited) LL that consists of 
two four-flux Pfaffians in each subband component. This 
state separates the usual monocomponent Pfafhan^’® in 
narrow quantum wells from the above-mentioned WC in 
wide quantum wells with low subband gaps. 


Our paper is organized as follows. We present the 
twocomponent model for wide quantum wells in Sec. 11. 
This model is used in the VMC analysis of the various 
quantum liquid phases in Sec. Ill, while Sec. IV is de¬ 
voted to the energy calculations of the different electronic 
crystals within the Hartree-Fock approximation. Because 
Secs. HI and IV are concerned solely with liquid and 
crystalline phases, respectively, they may be considered 
as intermediate steps for our main results for the overall 
phase diagrams in Secs. V and VI. They may therefore 
be skipped by readers interested by the resulting phase 
diagrams. However, they contain technical details and 
additional information that we consider sufficiently im¬ 
portant not to relegate them to an appendix. We discuss 
the theoretical phase diagram for the half-filled lowest LL 
in Sec. V in comparison with experimental hndings, and 
briefly discuss the phase diagram of the half-hlled second 
LL in Sec. VI. 


II. WIDE QUANTUM WELL MODEL 


Most experiments on quantum Hall effects are per¬ 
formed in Al^jGai-ajAs/GaAs quantum wells that are 
more than 40 nm The lowest subbands 

of the latter are very well approximated by the lowest 
subbands of an infinite square quantum well"^® For this 
reason we model the confinement potential by an infinite 
square quantum well of width w, in which the energy of 
the Qf-th subband is /{2mw‘^), where m is 

the effective mass of the electrons. The corresponding 
eigenfunctions are given by 

¥>o(C) = y^xJ “ (1) 

sin if a is even, 

where the ( variable is restricted to [—^, ^]. Although 
the Fermi energy usually lies below the first excited 
subband (a = 2) for non-interacting electrons, the N- 
particle ground state may be characterized by a non-zero 
occupation of higher subbands due to electron-electron 
interactions. Indeed, the presence of nodes in the excited 
subbands’ eigenfunctions lowers the interaction energy, 
and in very wide wells this energy gain may overcome 
the gap between the hrst and the second subband 
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in which case the excited subbands are populated. On the 
other hand, the subband energy quickly increases with 
the subband index a (Aq, oc a^), such that the ground 
state is unlikely to involve subbands beyond the first ex¬ 
cited one, at least for realistic well widths. In our model, 
we will only retain the lowest two subbands and justify 
this assumption a posteriori for the widest wells consid¬ 
ered. Hence we consider in the following a twocomponent 
system the two components of which represent the first 
{a = 1) and the second {a = 2) subband. Similarly to bi¬ 
layer quantum Hall systems, this twocomponent system 
may conveniently be described in terms of a pseudospin 
1/2, where the two eigenvalues of the pseudospin projec¬ 
tors correspond to the two subbands.In this analogy, 
the energy difference A between two subbands plays the 
role of an effective Zeeman effect that favors the occupa¬ 
tion of the lowest subband. Whereas one could in princi¬ 
ple also take into account the true spin of the electrons, 
within a four-component description, we rely on the fact 
that strong exchange interactions favor a full spin polar¬ 
ization at half-filling that is eventually oriented by the 
Zeeman field. The spin is thus effectively frozen out and 
omitted in the remainder of this paper. 

The three-dimensional Goulomb interaction 
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where r = (a;, y) is the position in the plane, gives rise to 
an effective twodimensional Coulomb interaction 





^ v^ai(C)y’ 0^2 (C')y^a3(C)y’a4(C') 


\/r2 + (C-C')^ 


(4) 


that depends on the four subband components ai, ..., 
04 through the integration over the associated wave func¬ 
tions (1). In the above expressions, Cq is the permittivity 
of the vacuum and e is the static dielectric constant of 
the host material. In the remainder of this paper, the 
energy unit is set to e^/dTreoe^s. 

When projected to the lowest LL, the full Hamilto¬ 
nian only consists of a Coulomb interaction term and a 
subband term. It reads 


H = V + — (ha=2 - na=i) (5) 

where the interaction operator is 

^ [ V'liWV'aaW 

X - r'\)'>pl^{r')ipaA^')(fr(fr' (6) 

in terms of the effective potential (4) and the twodimen¬ 
sional electronic field operators the lowest LL 

and subband a. Note that the usual normal ordering 
is taken into account by a proper regularization detailed 
in the appendix. Furthermore, ha = / 
is the number operator that counts the electrons in the 
subband a. 

For each FQHS, we compute the average value of the 
interaction energy (5) for a given number of electrons N 
by means of the variational Monte-Carlo method that 
is described in the appendix. The corresponding energy 
diverges in the thermodynamic limit, so that we choose 
to substract a regularization energy and divide by N in 
order to obtain a quantity that is finite even when ex¬ 
trapolated to IV —>■ oo; we call this quantity the cohesive 
energy. 

With regard to the subband degree of freedom, the trial 
states under scrutiny can be divided into two groups : 
monocomponent states, in which all electrons are in the 
same subband, and truly twocomponent states for which 
electrons necessarily occupy both components.Natu¬ 
rally, one expects monocomponent states (in a = 1) in 
the narrow-well limit, where the subband gap A is much 
larger than the interaction energy. The latter intervenes 
in the energy gain when populating the subband a = 2 
- because of the node in the wave function ipa=2{C)i the 
Coulomb repulsion is reduced in a = 2. From a theo¬ 
retical point of view, one may thus characterize a wide 
well as one with a width such that e^/dTrcoels is on the 
order of or larger than A. A more detailed discussion of 


the resulting scaling relations can be found in Sec. Ill C, 
where we analyze the phase diagram of the liquid phases. 

The notion of “components” does not necessarily re¬ 
fer to the subband index. Indeed, the components might 
refer to orthogonal combinations of the subband eigen¬ 
states. Here, we only consider the effective bilayer config¬ 
uration, in which the two components “-I-” and ” are 
defined by = {(pa=iiipa= 2 )/The effective bilayer 
wavefunctions are shown in Fig. 1. This choice for the two 
components is natural since the density profile of a cor¬ 
responding state has the z —z symmetry of the well. 
Our two-subband model thus accounts for the formation 
of an effective bilayer in wide wells due to the Coulomb re¬ 
pulsion, as corroborated by self-consistent Hartree-Fock 
calculations.^® However, we emphasize that our model is 
not biased to the formation of effective-bilayer states or 
the particular choice of ip± for the twocomponent basis. 
Indeed, the latter choice is natural only in the limit of a 
vanishing effective subband gap, where the renormaliza¬ 
tion of the gap by interactions is taken into account.^® 
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FIG. 1. (Color online) Lowest subband eigenfunctions (left) 
and effective bilayer wave functions (right). The latter are 
the -|- and — linear combinations of the former, i.e. = 
(tfii ± 1/P2) /%/ 2 . 


While monocomponent states may lie either in the low¬ 
est or in the first excited subband, all the twocomponent 
states that we will study are in the effective layer con¬ 
figuration, in which half of the electrons are in the </?+ 
state and the other half is in the (p- state. This as¬ 
sumption is corroborated by previous studies by Ref. 23 
that showed that the pseudospin-1/2 symmetry-breaking 
terms in the interaction Hamiltonian (3) favor an annihi¬ 
lation of the cc-component of the subband pseudospin, i.e. 
its orientation in the j/z-plane and a final orientation in 
the z-direction due to the subband gap. Furthermore, we 
have checked numerically in the present study that the 
lowest-energy monocomponent states have indeed a pseu¬ 
dospin polarization in the z-direction, whereas we have 
shown within exact-diagonalization studies via a varia¬ 
tional parameter that the twocomponent states occupy 
equally the p+ and p- states.^® 


III. LIQUIDS IN THE LOWEST LANDAU 
LEVEL: COMPOSITE-FERMION EERMI LIQUID, 
(331) HALPERIN AND PEAEEIAN STATES 


Before introducing the solid phases, we present our re¬ 
sults for the relevant liquid states. As mentioned above. 
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the monocomponent liquids, the CFFL and Pfafhan, can 
either be in the lowest or in the first excited subband. 

Concerning effective bilayer liquids, the (331) Halperin 
state is not the only one that needs to be considered. 
Indeed, beyond the (331) state we can also artificially 
superimpose monocomponent FQHSs, each with a filling 
factor 1/4, with a spinor structure. By doing so one can 
construct the “2CFFL” state, made out of two uncorre¬ 
lated CFFLs in each of the components, and the “2Pf” 
state built with two independent Pfafhan phases. These 
states are presented in detail in the following. 


A. Monocomponent liquids 


With the help of the variational Monte-Carlo (VMC) 
procedure described in the appendix A, we hrst compute 
the cohesive energies of the monocomponent liquids in 
the subband a. The trial wave function for the CFFL 2 n 


V'CFFLa^ {{a“}) = ^LLL det (^e 





i<j 


(7) 


can be expressed in terms of the Slater determinant 
det[exp(— -I- fc*z“*)], where Zj = x'^ + iy‘^ is the 
complex coordinate of the /-th electron in the subband 
a and kj the associated (complex) wave vector. Since 
the wave functions in the determinant bear non-analytic 
components, in terms of the conjugate coordinate z^*, 
the overall wave function needs to be projected into the 
lowest LL by the projector T^lll- Similarly, the Pfafhan 
reads 



where Pf denotes the Pfafhan of the matrix Mij = 
ll{zf — z^). In both wave functions (7) and (8), we have 

omitted the obiquous Gaussian factor exp(— kf P/4) 
for the lowest LL in order to simplify the notations. 
Moreover, they are both associated with the hlling factor 
iz = l/2n. 

For monocomponent states, we consider n = 1. The 
difference in cohesive energy of the CFFL 2 and the Pfaf¬ 
han is plotted in Fig. 2 as a function of the well width 
for both subbands. The Pf 2 state in the lowest subband 
{a = 1) is known to have a higher energy than the CFFL 2 
for ui = 0 in the lowest LL from previous studies. ^^43 pjg_ 
ure 2 shows that the CFFL 2 remains lower in energy than 
the Pf 2 at least up to = 20, although the energies 

get closer with increasing width. This piece of informa¬ 
tion is noteworthy since in wide quantum wells, the soft¬ 
ening of the Coulomb interaction leads to a reduced ratio 
V 1 /V 3 of the first two odd pseudopotentials compared to 
the pure Coulomb case, and this reduction is generally 


expected to be in favor of the Pfafhan state.As men¬ 
tioned above, the error bars for the energy difference are 
always larger than the uncertainty in the 1/N = 0 ex¬ 
traction shown in Fig. 15. 
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FIG. 2. (Color online) Cohesive-energy difference of the 
CFFL 2 and the Pf 2 in the lowest subband (a = 1) and in 
the first excited subband {a = 2) as a function of the well 
width wI£b- The CFFL 2 is always lower in energy. 


In the hrst excited subband (a = 2), although the 
energies of the Pf 2 and CFFL 2 are different from those 
computed in the lowest subband, their relative difference 
is almost the same (Fig. 2), such that also in the hrst 
excited subband the Pf 2 cannot be lower in energy than 
the CFFL 2 . One notices that the energies in the ex¬ 
cited subband are systematically lower than those in the 
lowest one (a = 1). This rehects the abovementioned 
node in the wavefunctions in the z-direction, which re¬ 
duces the Coulomb repulsion in a = 2 as compared to 
a = 1, and results in the tendency to spontaneously 
form bilayer states.^® However, we emphasize that one 
cannot compare states that belong to different subbands 
directly, since one has to specify the subband gap to 
do so. In Fig. 3, we have plotted the energy difference 
Acffl 2 = E^ffl 2 ~ -^CfIls between the CFFL 2 in the 
lowest and that in the first excited subband, neglecting 
the subband gap. The latter is plotted as an indepen¬ 
dent quantity for two values of the magnetic field (B = 5 
T and B = 10 T). Indeed, one notices that in the ab¬ 
sence of A, it is always favorable to form a CFFL in 
the first excited subband because of the lower interac¬ 
tion energy. If we now take into account the subband 
gap, the energy difference defines a gap value Acffl be¬ 
low which the CFFL prefers the first excited subband, in 
spite of the cost in interaction energy. In units of the typ¬ 
ical Coulomb energy e^/47reeo-^B the subband gap reads 
A ~ 5 . 93 VI 3 /{w/iB)^, such that the transition from the 
lowest subband to the first excited subband CFFL oc¬ 
curs at a width wjlB that depends on B. For instance 
for H = 5 T this transition occurs at w/is — 8 while it 
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happens at wjis ~ 9 for i? = 10 T (see Fig. 3). 
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FIG. 3. (Color online) Energy difference Acffl 2 = — 

EcFFh 2 between the CFFL 2 in the lowest and that in the first 
excited subband (red line). We have plotted the subband 
gap in the Coulomb energy unit Ab = A/(e^/47reeo^s) as a 
function of the width of the well w/£b for B = 10 T (green 
dashed line) and B = 5 T (blue dotted line) for comparison. 


B. Twocomponent liquids 

In addition to the above-mentioned monocomponent 
trial states, we investigate here the energies of two types 
of twocomponent states. The first one consists of states 
with no correlations between the components, namely 
the “uncorrelated CFFLs” ( 2 CFFL 4 ), which consist of 
two CFFLs at = 1/4 in each component, and the 
“two uncorrelated Pfafhans” ( 2 Pf 4 ) composed of four- 
flux composite fermions. In terms of the complex posi¬ 
tions of the /-th electron in the effective-bilayer state 
±, the 2 CFFL 4 wave function reads 

'02CFFL4 ^■2^1 ) ■ • ■ ) J ' ■ ' > ^N/2 

= i’CFFhi > ■ ■ • ) tpCFFLi [z^ , • • ■ , %/2) 

where ' 0 cffl 4 , given by Eq. (7) for n = 2, is the four-flux 
CFFL wave function. The 2 Pf 4 wave function is defined 
analogously as a product of two four-flux Pfafhan states 

V'2Pf4 ({^i }) ~ '4’Ffi (^Z^, ■ ■ ■ , z'1^^2^ '4’P{i (^4 , . . . , Zj^^^ , 

( 10 ) 

where is defined by Eq. ( 8 ). Finally we investigate 
the (331) Halperin state^° 

^(331) [zt,..., z"/2) 

=n ~ 4 f n n (4 - z~) 

i<j k<l m,n 


in which the electrons in both components are now cross- 
correlated by the last term. Similarly to the monocom¬ 
ponent wave functions (7) and ( 8 ), we have omitted in all 
trial wave functions (9), (10) and (11) the Gaussian factor 
exp[-(X]f’^ \zt\‘^ + J2f~ \zi?)m for the lowest LL. No¬ 
tice that one can also construct another Halperin state, 
(113) with inverted exponents, at p = 1/2. This state 
is however unstable and subject to phase separation^^ so 
that we do not consider it here. 

The cohesive-energy differences of those states are 
compared in Fig. 4. We have chosen the cohesive en¬ 
ergy of the (331) state as the reference. The latter has 
the lowest energy for w/£b < 15 while the 2 CFFL 4 is fa¬ 
vored in very wide wells (wjiB > 17). The precise value 
of the width at which the transition occurs cannot be ex¬ 
tracted from our numerical results due to uncertainties 
associated with the VMC and the polynomial extrapola¬ 
tion (see inset of Fig. 4). This result can be understood 
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FIG. 4. (Color online) Cohesive energies of the effective-layer 
liquid states. The reference chosen here is the cohesive energy 
of the (331) state. The inset is a zoom around the transition 
between the (331) state and the 2 CFFL 4 . Notice that the 
error bars prevent the precise location of the transition (be¬ 
tween w/£b — 15 and w/Ib — 17). 


quite intuitively: in a wide quantum well, correlations 
between the different components become less relevant, 
similarly to bilayer systems in the large-distance limit. 
One therefore retrieves two independant onecomponent 
systems with = 1/4, where the ground state is the 
2 CFFL 4 ."‘®d 7 jygi; as in the monocomponent case the 
2 Pf 4 state is systematically higher in energy than the 
2 CFFL 4 and the (331) state, hence according to our cal¬ 
culations no Pfafhan phase can be stabilized in the lowest 
LL in a wide quantum well. 
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C. Phase diagram of mono- and twocomponent 
liquid states 


To be relevant for realistic systems, the phase di¬ 
agram of the liquid states should take into account 
the subband gap A. Indeed, it adds a contribution 
^™bband = NA/2 to the effective bilayer states and a 
contribution if|ubj 3 and = the excited subband 

states. The phase diagram is shown in Fig. 5 as a func¬ 
tion of well width and subband gap. The subband gap is 
treated here as a parameter independant of the width for 
illustration. As expected, large subband gaps favor the 



FIG. 5. (Color online) Phase diagram for liquid states. Black 
indicate regions where the uncertainties of the lowest-energy 
state overlaps with the uncertainty of another state, such that 
we cannot conclude. No Pfafhan state is stabilized here. Note 
that the black region of undeterminacy that separates the 
2 CFFL 4 and (331) phases is rather large since the two phases 
have close energies for w/£b — 15, and they have the same 
subband energy per particle F,^ubband = A/2. 


lowest subband CFFL 2 since the latter is the state with 
lowest energy in monocomponent systems. Conversely, 
when the well is large and for moderate subband gaps, it 
becomes favorable to populate the first excited subband 
due to the gain in interaction energy. In that regime 
(lower right part of the diagram), the CFFL 2 is also 
stabilized but now in the first excited subband. These 
two regimes of monocomponent states are separated by 
a region in which the subband gap is on the same or¬ 
der of magnitude as the gain in interaction energy when 
populating the first excited subband, such that an effec¬ 
tive subband gap that one can define on the mean-field 
level vanishes. In this approach, one finds for this re¬ 
gion the scaling behavior for the correlation energy gain 
A 2 -comp of the twocomponent states compared to the 
lowest subband^^ 


A2—comp 

e^/dTreoefs £b 


( 12 ) 


that is in good agreement with our results depicted in Fig. 
(5). In this region, one stabilizes the various twocompo¬ 
nent states of the effective bilayer, i.e. the 2 CFFL 4 if the 
well is very large {w/£b > 17) and the incompressible 
(331) state for w/£b < 15. 


IV. WIGNER CRYSTALS IN WIDE QUANTUM 
WELLS IN THE LOWEST LANDAU LEVEL 

The FQHSs presented in the previous section are in 
close competition with differents types of Wigner crys¬ 
tals (WCs).^®’"^® In contrast to the FQHSs, the latter 
are charge density waves of a particular type that can 
be described conveniently within the Hartree-Fock ap¬ 
proximation of the interaction Hamiltonian.They 
are relevant since, for instance, the triangular WC has a 
lower energy than all known FQHS at low densities'*®’^^ 
V < 1/6, and it also shows up in the vicinity of FQHE 
magic filling factors through the reintrant integer quan¬ 
tum Hall effect.^®’®® 


A. Monocomponent Wigner crystals 

Let us first concentrate on monocomponent WCs in 
one of the subbands a. In WCs the electrons minimize 
the Coulomb interaction by being situated at the nodes 
of a lattice. The localized wave function of an electron 
around R with lowest kinetic energy is a coherent state 

d'R(r) = (r|e^^V,0) = 

V 27r 

(13) 

and the corresponding many-body state involves anti- 
symetrization A of N localized states, where N is the 
number of electrons 


'kwc (ri, • ■ • An) = A ['kRi(ri)... 4'R„(rAr)] . (14) 


If we neglect the overlap of neighboring localized states, 
the probability density of the many-body state (14) is 


p(ri) = J dV2...d^rAr|4'wc(ri,---,i'Af)l^ 



j 




(15) 

where p is an element of the permutation group ctn- Ne¬ 
glecting the overlap between adjacent lattice sites is a 
reasonable assumption in the low-density limit, where 
the average distance between electrons is much larger 
than the typical extent £b of the wavefunctions. In the 
present case, however, we investigate WCs at v = 1/2, 
where overlap corrections become relevant. In this case. 
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the square of the probability density in reciprocal space 
that is more convenient for our Hartree-Fock calculations 
reads 

|(p(q))P= e-^^/^|(p(q))P, (16) 

whose explicit expression is derived in Appendix B [Eq. 
(B27)]. 

In terms of the modulus square of the projected den¬ 
sity I (p(q)) P, which plays the role of the crystalline or¬ 
der parameter, the energy of WCs in the Hartree-Fock 
approximation reads^^ 

^HF = l(p(q))l^ (17) 

q 

where A is the area of the system, and ■!;hf(9) = '5h( 9) — 
■ijF(Q) is the Hartree-Fock potential. The direct (Hartree) 
term 

VH(q) = v(q) exp(-qV2) (18) 

is the lowest-LL-projected Fourier transform of the po¬ 
tential V{r), and the exchange (Fock) term 

VF^q) = vh {q£%) (19) 

is proportional to the (inverse) Fourier transform of vh 
for a rotationally invariant interaction potential. Notice 
that the formalism is readily generalized to higher LLs 
with index n (n = 0 being the lowest LL) via the re¬ 
placement of the Gaussian exp(—5^/2) in Eq. (18) by the 
complete form factor exp(—g^/2)[L„((jf^/2)]^, in terms of 
the n-th Laguerre polynomial. 

The density riei = N/A is related to the filling factor iz 
by riel = ly/2Tr£'^, so the energy per particle finally reads 

^ = ( 20 ) 

q 

In the lowest subband the relevant potential is [see 

Eq. (4)], it is I /2222 gj.g|. g^cited subband, and in 

the effective bilayer configuration we have to specify both 
the intra-layer interaction and the inter-layer one. As 
mentionned in Sec. HI (and discussed in the appendix) 
the energy of Eq. (20) needs to be regularized to ob¬ 
tain the cohesive energy Ecoh.- We have used the same 
procedure as for the one for the liquid state (see the ap¬ 
pendix). The g = 0 term of the sum in the energy (20) 
is (partly) canceled by the regularization energy, and the 
rest of the sum converges rapidly thanks to the Gaussian 
factors exp(—g^/2), such that only a hundreds of terms 
are enough to reach machine precision. 

We now have to determine the lattice that minimizes 
the energy of the crystal for both of the lowest subbands. 
We consider a general (monoclinic) twodimensional re¬ 
ciprocal lattice and minimize the cohesive energy with 
respect to the vectors that define the elementary cell, 
while keeping the area of the unit cell fixed to I/nei. In 
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FIG. 6 . (Color online) Anisotropy parameter a of the rect¬ 
angular Wigner crystal with lowest energy, in the lowest and 
first excited subband, as a function of the well width w. The 
lattice becomes more anisotropic with increasing width. 

doing so we find that the minimum is obtained for a rect¬ 
angular lattice, in both subbands and for all well widths. 
The basis vectors of this rectangular lattice read 

al = ^={^/a,0) and a 2 = —^(0,1/Va), 
^/nel WAi 

( 21 ) 

where we have introduced the anisotropy parameter a > 
1 such that jjaill = Q!||a 2 ||. With this choice, the lat¬ 
tice sites are approached in the y-direction, whereas 
the lattice spacing is increased in the ^-direction. The 
anisotropy parameters of the crystal with lowest cohesive 
energies are shown in Fig. 6. The Wigner crystal with 
lowest energy appears to have a non-negligible anisotropy 
even in thin wells, and this anisotropy grows with in¬ 
creasing well width. Anisotropic WGs are peculiar to 
half-filled LLs, and were interpreted as cousins of stripe 
phases in an earlier work (in the third LL).^® Interest¬ 
ingly, we find even in the lowest LL (if we artificially dis¬ 
card the quantum-liquid phases) a tendency of the WG to 
form stripe-like phases. Notice that we have considered 
also the usual triangular WG, but it is higher in energy 
than the rectangular one. 


B. Twocomponent Wigner crystals 

As for the liquid phases, in addition to the monocom¬ 
ponent states that are either entirely contained in the 
lowest or in the first excited subband, we also have to 
consider effective bilayer WCs that consist of a WG in 
each of the effective layers. In order to minimize the in¬ 
terlayer repulsion the WG of one layer is positioned at 
the interstitials of that in the other layer. Here, we only 
present two types of lattices for the two component WG, 
the stacked square and stacked triangular lattices that 








(a) Stacked square 



(b) Stacked triangular 
,Q.-P 



FIG. 7. The stacked square and stacked triangular lattices. 
The lattice of the — layer (empty circles) is shifted with re¬ 
spect to the one of the -I- layer (filled circles) by a displacement 
vector a. 


are depicted on Fig. 7. In those WCs half of the elec¬ 
trons are in the -I- effective layer, while the other half is 
in the — layer and localized around the sites of a lattice 
that is shifted by a displacement vector a with respect 
to the one of the + layer. We have also considered more 
general stacked WCs such as the stacked rectangular one, 
but they happen to be higher in energy contrarily to the 
monocomponent crystals. 

The intra- and inter-layer Hartree potentials are given 
by 


4^/^)(q) = i}(A/E)(g)e-4 (22) 

where and are the Fourier transforms of the 
intra-layer effective potential (A5) and inter-layer effec¬ 
tive potential (A6), respectively. Careful handling of 
the Hartree-Fock approximation shows that the corre¬ 
sponding Fock terms are v^\q) = 2 tt£b v^\q£g) and 
vl^\q£'^), such that the inter-layer Fock 

/•p\ 

potential Vp is related to the cross term (A7) and not 
to the inter-layer potential (A6). 

In terms of those potentials the Hartree-Fock energy 
per particle of a stacked WC reads 



0 5 10 15 20 25 

wI£b 


FIG. 8. (Color online) Cohesive energy difference of 

the stacked triangular and stacked square Wigner crystals as 
a function of the well width w. The stacked square crystal has 
lower energy for wjiB < 13, while the stacked square lattice 
is favored in larger wells. 


displacement of the two triangular lattices in each com¬ 
ponent. Conversely, when the effective layers get closer 
to each other, the lattice configuration tries to accomo¬ 
date simultaneously the inter- and inter-layer interac¬ 
tions. Therefore the overall lattice, which can be viewed 
as a bipartite lattice with one sublattice per layer, should 
be as closely packed as possible. This is not the case of 
the two interpenetrating triangular lattices that form a 
honeycomb lattice, and a stacked square lattice is there¬ 
fore favored. Notice that an overall triangular lattice, 
which would have the closest packing, cannot be formed 
at finite separations between the effective layers because 
it is not bipartite. It might eventually be realized in the 
putative limit where the two effective layers are no longer 
spatially separated and where the intra- and interlayer 
repulsions are equal. 


C. Phase diagram of mono- and twocomponent 
Wigner crystals 


A. 


stacked 

Iv 


27r ^ 
q 


( 9 ) + 41^^ ( 9 ) cos(q • a) I (p(q)) 1^ 


(23) 

where the reciprocal lattice that defines (p(q)) accomo¬ 
dates only half of the electrons, and a is the relative dis¬ 
placement vector of the stacked lattices (see Fig. 7). 

The energies of the two crystals are compared in Fig. 8, 
where we plot the difference in cohesive energy between 
the stacked triangular and the stacked square crystals. 
The results can be understood easily in comparison with 
the true bilayer quantum Hall system. Similarly to the 
latter,^® the stacked triangular crystal is favored for large 
distances between the effective layers {i.e. large well 
widths). Indeed, when the effective layers are well sepa¬ 
rated the lattice tries above all to minimize the intra-layer 
repulsion, which favors triangular lattices. The residual 
inter-layer repulsion is then accomodated by the relative 


The WCs phase diagram presented in Fig. 9 shares 
common features with the one of liquids (Fig. 5). In 
large wells and for small subband gaps A the excited 
subband is occupied due to the reduced interaction en¬ 
ergy, whereas in thin wells and for large subband gaps, 
a rectangular WC in the lowest subband has the low¬ 
est energy. Between these two limits of monocomponent 
states, one finds a region where stacked WCs can be sta¬ 
bilized (notably only in wells wider than 1 Q£b), in the 
square version for wj £b ^13, and in the triangular one in 
the top right orange region of the diagram. The stacked 
WCs are indeed favored by inter-component correlations 
and are thus reminiscent of the (331) liquid state. As 
for the twocomponent states in the liquid phase diagram 
(Fig. 5), the stacked crystals are encountered in a region 
that roughly obeys the scaling behavior shown in Eq. 
( 12 ). 
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FIG. 9. (Color online) Wigner crystals phase diagram. The 
parameters are the width of the well w and the subband gap A 
(here in units of e^/47reeo^s)- “Rec.” denotes the rectangular 
WC in the lowest subband, “Rec. exc. sub.” is its alter ego 
in the first excited subband, “St. Tr.” designs the stacked 
triangular WC and “St. Sq.” stands for stacked square. 


V. PHASE DIAGRAM OF THE HALF-FILLED 

LOWEST LANDAU LEVEL - COMPARISON 
BETWEEN LIQUID AND SOLID PHASES 

In the past two sections, we have presented an ex¬ 
tensive set of model states for the 2DES at v = 1/2. 
Whereas the energies of the various liquid states were 
computed by means of the VMC method in the spherical 
geometry, the energies of the solid states were obtained in 
the Hartree-Fock approximation in the planar geometry. 
We now give a detailed comparison of all investigated 
phases and draw a phase diagram for the competing liq¬ 
uid and solid states. The spherical and disc geometries 
are expected to lead to the same energies in the thermo¬ 
dynamic limit, which can thus be compared directly. The 
resulting phase diagram is shown in Fig. 10 as a function 
of the well width w/is and the subband gap A, which we 
continue to keep as an independent parameter for the mo¬ 
ment. Whereas for narrow wells and large subband gaps, 
we retrieve, as expected from previous studies,liq¬ 
uid ground states in the half-hlled lowest LL, these liquid 
phases disappear in favor of a rectangular WC in the first 
excited subband in the limit of large wells and small sub¬ 
band gaps. In addition to the usual CFFL 2 state in the 
lowest subband, the (331) Halperin state remains stable 
in a part of the intermediate region, where both sub¬ 
bands are occupied even if it occupies a less prominent 
part of the phase diagram as compared to that (Fig. 5) 
where only liquid phases were considered. Note that the 
2 CFFL 4 state may also be stabilized in this intermediate 
region, albeit for A > 0.3 and w/ts > 13, i.e. for narrow 
wells and very large magnetic fields that are not realized 
experimentally. 

In order to make a closer connection with experiments. 


FIG. 10. (Color online) Phase diagram of the wide quantum 
well a.t u — 1/2. Here the parameters are the width w in units 
of the magnetic length Ib and the subband gap A in units of 
the typical Coulomb energy ^ j^-Kteal b- 

we need to plot the phase diagram in Fig. 10 in terms 
of experimentally relevant parameters. In contrast to 
our theoretical description, the subband gap in units of 
e^/ 47 reeo^B and the well width in units of the magnetic 
length cannot be varied independently. First, the well 
width is fixed by the sample-fabrication process, and one 
therefore needs to modify the ratio ic/fs oc by vary¬ 
ing the magnetic field B, i.e. an increase of the magnetic 
field drives the system towards the large-well limit. More¬ 
over, in order to keep the filling factor v = n^xj(eBjK) 
constant when varying S, one has to change the elec¬ 
tronic density riei simultaneously by applying an external 
gate voltage to the sample. An increase in B, however, 
also raises the energy unit e^/AneeoiB oc VB, thus reduc¬ 
ing effectively the subband gap A in our units. Hence for 
a given width w, increasing B (and consequently riei) re¬ 
sults in a larger width w/Ib and a reduced subband gap 
A in our phase diagram (Fig. 10). 

The phase diagram as a function of realistic values for 
the well width w and the electronic density riei is plotted 
in Fig. 11, in comparison with a phase diagram by Sha- 
bani et for measured compressible, incompressible 
and insulating phases (Fig. 12). As mentioned above, at 
a fixed value of the well width, the low-density limit cor¬ 
responds to rather low magnetic fields where the subband 
gap A is substantial when compared to the interaction 
energy scale. The limit of large subband gaps and small 
well widths is therefore located in the lower left corner of 
the phase diagram, whereas the large-ru/small-A limit is 
found in the upper right part. The scaling behavior (12) 
of the intermediate region separating these to limits now 
translates to w oc l/riei oc l/H in the phase diagram in 
Fig. 11. Furthermore, as mentioned above, one notices 
that the only relevant twocomponent phase in this inter¬ 
mediate region is the incompressible (331) Halperin state 
the 2 CFFL 4 being shifted to parameters that are not in¬ 
vestigated experimentally. Before comparing our findings 
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FIG. 11. (Color online) Theoretical phase diagram as a func¬ 
tion of the experimentally relevant parameters, the width w 
in nm and density rigi in units of 10^^ cm“^. 


in more detail with the experimental phase diagram (Fig. 
12 ), we point out that we have checked explicitly, within 
an enlarged model taking into account the third subband, 
that the latter is only populated for very large wells that 
are usually not studied experimentally {w > 100 nm or 
riel ^ 4.5 X 10^^ cm“^), i.e. beyond the parameters for 
which we have presented the phase diagrams. 

In experiments one does not have direct access to the 
nature of the ground state, but only to its transport char¬ 
acteristics. The experimental phase diagram in Fig. 12 
from Ref. 21 shows that a compressible state in the limit 
of narrow quantum wells and low electronic densities is 
separated from an insulating state at large w and high 
densities via an intermediated region where an incom¬ 
pressible FQHE is observed. This is in excellent agree¬ 
ment with our phase diagram, where the FQHE is at¬ 
tributed to a (331) state and the compressible one to a 
CFFL 2 . Notice that even the lines separating the differ¬ 
ent regions are in good quantitative agreement. Whereas 
the compressible and incompressible states have com¬ 
monly been understood in terms of a CFFL 2 and the 
(331) state, respectively,^^’^^ the origin of the insulating 
state remained unexplaned even if a WC seemed a natu¬ 
ral candidate in view of the pertinence of correlations at 
half filling. Our results indicate that the measured insu¬ 
lating state can indeed be understood in terms of a crys¬ 
talline state. Most saliently, it is a rectangular WC that 
is formed in the first excited electronic subband due to 
the confinement potential. Our results indicate that the 
effective interaction potential is modified by the subband 
wave functions in such a manner that crystalline rather 
than liquid phases are favored there, whereas crystalline 
phases at half filling are usually encountered in higher 
LLs in the form of stripe phases.Notice furthermore 
that the rectangular shape of the WC in the first excited 
subband is reminiscent of the stripe phase when consid¬ 
ering the electronic density, especially since the latter in 


higher LLs reveals an unstability towards a density mod¬ 
ulation along the stripes. 


B,.,RT) 



FIG. 12. (Color online) Experimental phase diagram of the 
wide quantum well at z/ = 1/2, taken from Ref. 21. Here W 
in the width of the well, and n is the electronic density. 


VI. PHASE DIAGRAM OF THE HALF-FILLED 
SECOND LANDAU LEVEL 

In this last section, we present the phase diagram for 
electron solids and quantum liquids for the half-filled sec¬ 
ond LL. From a technical point of view, the difference 
between the two LLs arises from the difference in the ef¬ 
fective interaction potentials, which take into account the 
orbital wave functions associated with the LL. Concern¬ 
ing the VMC calculation of the liquid phases in the sec¬ 
ond LL, we use a different interaction than in the lowest 
one. Whereas, in the lowest LL, we could use a polyno¬ 
mial interaction (see the appendix) that reproduces ex¬ 
actly the pseudopotentials of the effective interaction (4), 
we have no simple formula to compute the pseudopoten¬ 
tials in the second LL associated with monomials on 
the sphere. Thus, similarly to Ref. 43, we use an alter¬ 
native Gaussian interaction. The coefficients are again 
chosen such that the associated pseudopotentials, now 
computed in the planar geometry, match the first five 
odd pseudopotentials of the effective interaction (4). 

Figure 13 shows the energy differences for some of the 
states, with respect to the (331) state. Here, we have 
chosen, for illustration reasons, to show only intrinsic 
twocomponent states that are insensitive to the subband 
gap. One notices that the stacked WC remains higher 
in energy than the different twocomponent FQHSs. Fur¬ 
thermore, we have checked that a phase of alternating 
stripes in the effective-bilayer configuration is even higher 
in energy than the stacked WC. Therefore in the interme¬ 
diate region that separates the different monocomponent 
phases, similarly to the phases in the lowest LL, one finds 
a competition between the (331) state, which is favored 
in narrow wells, and the 2 Pf 4 in wide quantum wells, 
with a transition at wjlB ~ 3. We emphasize that the 





11 



FIG. 13. (Color online) Energy differences of various in¬ 
trinsic twocomponent states, as a function of the well width. 
We have chosen the (331) state as the reference state (black 
continuous line). 

incompressible 2 Pf 4 state is a specificity of the second LL 
and is not stabilized in the lowest one. 



g 
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FIG. 14. (Color online) Phase diagram of the various states 
in the second LL as a function of the electronic density Uei 
and the well width w. The (Pf 2 ) encountered phases are the 
monocomponent Pfafffan in the narrow-well limit, the rect¬ 
angular WC in the limit of wide quantum wells, whereas an 
incompressible state (2Pf4) of two four-flux Pfaffians (one in 
each subband component -|- and — is stabilized in an inter¬ 
mediate regime.) 

The phase diagram for the phases in the second LL is 
depicted in Fig. 14, as a function of the experimentally 
relevant parameters w and riei- One notices that of the 
above-mentioned twocomponent states in the intermedi¬ 
ate region of vanishing effective subband gap only the 
incompressible 2 Pf 4 phase survives, whereas the (331) 
state is covered by the usual monocomponent Pfafhan 
Pf 2 in the narrow-well limit, which extensive theoret¬ 
ical work^^^^^’"^^ showed to be at the origin of the in¬ 
compressible FQHS observed in this limit.Notice fur¬ 
thermore that the Pfafhan origin of the 5/2 FQHS has 
recently been discussed intensely based on effective in¬ 


teraction potentials modihed by LL mixing that breaks 
the particle-hole symmetry of the half-hlled LL.®^“^^ As 
a consequence, it was argued that the anti-Pfafhan (the 
hole version of the Pfafhan) would be a more appropriate 
groundstate wavefunction at this hlling factor. 

However, a possible transition between Pfafhan and anti- 
Pfafhan sensitively depends on the degree of LL mixing, 
eventually altered by hnite width effects, and for some 
of these parameters the Pfafhan is indeed stabilized. 
While this debate is still ongoing, we emphasize that we 
neglect LL mixing within our model, such that the mono¬ 
component states are particle-hole-symmetric, and Pfaf¬ 
han and anti-Pfafhan therefore have the same energy. 
Whereas the FQHSs are therefore different from those in 
the lowest LL, we obtain the same rectangular WC in 
the hrst excited subband as the ground state in the limit 
of wide quantum wells and large densities, in agreement 
with the insulating state in the lowest LL. However, we 
emphasize that this WC is encountered at larger values 
of w and riei than in the lowest LL. 


VII. CONCLUSION 

In conclusion, we have compared the energies of various 
model states for the 2DES at z/ = 1/2 in the lowest LL, 
taking into account the hnite width of the conhning quan¬ 
tum well explicitly within a square quantum well model. 
Keeping only the lowest two subbands, we obtain an ef¬ 
fective twocomponent model of electrons restricted to a 
single LL in which the subband wave functions modify 
the effective interaction potential. This approximation 
is justihed for the investigated parameter range, which 
covers that of available experimental data and where we 
have checked explicitly that higher subbands beyond the 
second one are unoccupied. Within our twocomponent 
model, we have investigated in a variational approach 
the energies of various candidates for the ground state 
when varying the well width as well as the subband gap. 
This approach consists of VMC calculations for the en¬ 
ergies of the different liquid states, such as the compress¬ 
ible CFFL, twocomponent CFFL ( 2 CFFL 4 ) and the in¬ 
compressible (331)-Halperin, Pfafhan and twocomponent 
Pfafhan ( 2 Pf 4 ) states, that we compare to Hartree-Fock 
calculations for the energies of the electron solids, namely 
the various mono- and twocomponent WCs. 

Globally, one hnds two limiting regions that can be 
seen as effective monocomponent limits. The hrst one is 
that of low well width and large subband gaps in which 
case all electrons reside in the lowest electronic subband. 
In this limit, which has largely been studied in the liter¬ 
ature, one retrieves the compressible CFFL in the lowest 
LL. The second (monocomponent) limit is that of large 
well widths and low electronic subband gaps, where it is 
energetically favorable to populate the excited subband. 
The energy cost due to the subband gap is then over¬ 
come by the gain in interaction energy since the wave 
functions of this subband have a node that effectively re- 
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duces the Coulomb repulsion as compared to the lowest 
subband. Somewhat unexpectedly, one finds that crys¬ 
talline phases are favored over the usual quantum liquids 
in this limit. Most saliently, we obtain a rectangular WC 
in the excited subband as the ground state in this limit. 
This phase is consistent with the experimentally observed 
insulating electronic behavior in this limit. 

The above-mentioned limiting regions of monocompo¬ 
nent states are separated by a relatively narrow region 
in which the subband gap is on the same order as the 
interaction-energy gain in populating the excited sub¬ 
band. In this region the equal population of both sub¬ 
band states, or more precisely of the symmetric and anti¬ 
symmetric bilayer-type superposition of both, allows for 
a reduction of the correlation energies. One hnds then 
intrinsic twocomponent states, both if one considers crys¬ 
talline phases, in the form of stacked crystals, and liquid 
phases, in the form of the (331) Halperin state and a state 
consisting of two uncorrelated four-flux CFFL ( 2 CFFL 4 ) 
in each of the two components. A comparison between 
these different phases shows that the liquid states are fa¬ 
vored over the stacked crystals and that the incompress¬ 
ible (331) state is realized in the region of experimentally 
accessible parameters. Our theoretical phase diagram 
is in excellent agreement with experimental findings for 
the different compressible, incompressible and insulating 
phases. First, our study confirms that one can under¬ 
stand the compressible phase in terms of the CFFL and 
the incompressible one as a (331) state. Second it identi¬ 
fies a rectangular WC in the hrst excited subband as the 
origin of the insulating phase in the limit of large well 
width and low subband gap. From an experimental point 
of view, it is challenging to verify our image of a rectangu¬ 
lar WC in the first excited subband as compared to other 
possible (crystalline) phases. In principle, a measure of 
the charge density in the ^-direction would be required 
that is beyond the scope of present experimental tech¬ 
niques. Notice that recent microwave experiments have 
been interpreted in terms of a crystalline phase such as 
ours,^^ but an identification in terms of the particular 
crystal we obtain would require a detailed study of the 
pinning mode, which might bear information about the 
charge profile in the z-direction. Furthermore, one may 
speculate that at intermediate temperatures the rectan¬ 
gular WC partially melts into a stripe-type phase restor¬ 
ing partially translation symmetry along the stripe direc¬ 
tion. This would manifest itself in a highly anisotropic 
magnetoresistance as such measured in higher LLs.^°“^^ 
However, these possible effects are mere speculation at 
this stage, and a detailed study of both the pinning mode 
and a possible partial melting is beyond the scope of the 
present paper. 

Finally, we have investigated the competition between 
electron liqnid and solid phases in the half-filled sec¬ 
ond LL. Similarly to the lowest one, monocomponent 
states in the narrow well are separated from monocompo¬ 
nent states in the wide-well small-subband-gap limit by 
an intermediate region of intrinsic twocomponent states. 


However, whereas one obtains the same insulating phase, 
which originates from a rectangular WC in the first ex¬ 
cited subband, at large values of w and riei, the FQHSs 
are different from those in the lowest LL. Instead of 
a competition between a monocomponent CFFL and a 
twocomponent (331) Halperin state, we encounter a com¬ 
petition between two Pfafhan phases in the second LL. In 
the narrow-well limit, we are confronted with the usual 
monocomponent Pfaffian state, which vanishes in the in¬ 
termediate region in favor of the 2 Pf 4 state that consists 
of two four-flux Pfaffian states in each effective-bilayer 
component. Whereas both liquid phases are incompress¬ 
ible, the WC in wide wells is expected to manifest itself 
via an insulating phase. 
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Appendix A: Variational Monte-Carlo 


In this appendix we present the procedure we have 
used to determine the energies of the various liquid model 
states considered in this paper. 

We compute the energies of the relevant FQHSs by 
direct integration of the trial wave functions in real space. 
For a system of N electrons, the interaction energy in a 
model state 'k reads 


{V)q, =-^ J ... (fvN H (ri,..., Tat) 


X |5'(ri,...,rAr)| , (Al) 


where 
Af = 


J (fn ... (frN (ri,..., r 7 v)|" 


1 -1 


(A2) 


is a normalization constant that is necessary since we deal 
with unnormalized wave functions. The integrals (Al) 
and (A2) are computed by VMC in the spherical geom¬ 
etry. 

Since we do not use a normal-ordered version of the in¬ 
teraction operator in Eq. ( 6 ), the interaction energy per 
particle diverges in the thermodynamic limit. We there¬ 
fore have to introduce an explicit regularization energy 
that is common to all trial states in order to allow for a 
direct comparison between the latter. We arbitrarily set 
this regularization energy Areg. to the homogeneous part 
of the interaction energy of an electron system that only 
populates the lowest subband 

Areg. = ^ y d'^r. (A3) 
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This energy also diverges in the thermodynamic limit, 
but if we substract it from the interaction energy (Al) the 
interaction energy per particle becomes finite. We call it 
the cohesive energy E^oh., and will use this quantity in 
the following to compare the different phases, be they 
liquid or solid. We have checked that other choices for 
the reference energy do not change the energy differences 
between trial states. 



l/N 


FIG. 15. (Color online) Example of extrapolation of the 
cohesive energy to the thermodynamic limit. Here the state is 
the CFFL in the lowest subband, and the width is w = 10 
Each blue point was obtained by averaging 10 VMC runs of 
10 ® iterations, hence the smallness of the error bars, which 
is of the size of the blue symbols. Since we consider only 
filled shells, we only access values of N that are square of an 
integer. The solid red line is the result of the linear regression, 
and the dashed green line shows the second order polynomial 
regression. The cohesive energy is the value of the second- 
order polynomial for l/N = 0, and the error is the difference 
of this value with the one given by the linear regression. 


For all accessible sizes we first compute the interaction 
energy (Al) by VMC, then we substract the regulariza¬ 
tion energy (A3), which is finite for a hnite size of the 
sphere. A division by the number of electrons N gives 
the finite-size cohesive energy. Then we plot this quan¬ 
tity as a function of l/N and perform a second-order 
polynomial regression. The cohesive energy in the ther¬ 
modynamic limit is given by the value of this polynomial 
at l/N = 0. The procedure is illustrated in Fig. 15, 
where a comparison with a linear regression to l/N = 0 
shows that the uncertainty (< 0.1 %) of the procedure is 
sufficiently small with respect to the energy difference of 
the trial states. 

In the VMC method the interaction potential (4) is 
evaluated many times. In order to speed up computa¬ 
tions we have used an alternative potential defined by 


14>oiy. (^) 


1 

h 



(A4) 


is the radius of the sphere for the system size under con¬ 
sideration. Here, the coefhcients Cfc are chosen to be such 
that the Haldane pseudopotentials of the polynomial in¬ 
teraction (A4) reproduce exactly the ones of the original 
interaction (see Ref. 23 for details). Since the pseupoten- 
tials are identical, so are the interaction energies, such 
that this method does not change the computed energy 
values. 

In order to analyze the twocomponents states, we have 
to determine two sets of coefficients Cfc for the polynomial 
interaction coefficients of Eq. (A4) defined above, one for 
the intra-layer interaction 


y(A)(^) 



|yp+(C)ri^+(OI^ 

+ (C - CO^ ^ 


(A5) 


between electrons that belong to the same effective layer 
[the effective intra-layer interaction is the same in both 
layers due to <p+(C) = and another set for the 

inter-layer interaction 


V(®)(r) 




+ (C-CO^ 


(A6) 


In principle one should also take into account the “cross 
term” 


W W 

v(^)(r)= r r ^+(C) ^-(C0y^-(C)^+ (O^ 

i-f i-f + (C - C'Y 

(A7) 


This term vanishes for a real bilayer but not for an ef¬ 
fective bilayer formed in a quantum well, due to the non 
zero overlap of the -I- and — effective layer states. How¬ 
ever, 1/1^1 is more than an order of magnitude smaller 
than or V® (even for the largest well widths con¬ 
sidered here), such that we neglect it in our variational 
calculation. 

For completeness and a check of our numerical meth¬ 
ods, we present in Tab. I our results for the liquid phases 
in the limit ic = 0 of the ideal 2DES in comparison with 
published results from Refs. 43 and 27. The results agree 
within the numerical error bars. 


State 

Ref. 43 

Ref. 27 

This work 

(331) 

— 

-0.4631(3) 

-0.4633 ±0.001 

Pfafhan 

-0.4569(2) 

-0.4573(3) 

-0.4574 ±0.003 

CFFL 

-0.46557(6) 

— 

-0.4656 ± 0.004 


TABLE I. Cohesive energies of the liquid states we consider 
in the w = 0 “pure Coulomb” limit. The first two columns 
were retrieved from existing publications, and the last column 
gives our results in this limit. 


where is the number of flux quanta and R 
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Appendix B: Overlap corrections in the Wigner 
crystals 

1. Preliminary: 2-body localized state 


and 


Sp{r) = 


AA 2 

- cos , „ 

TT V 2 


/ yRp 




(B9) 


A normalized coherent state centered around R reads 


4'r (r) = (r|R) = 




e 2 


iiyX-xY) 


(Bl) 


Before writing the state for two electrons localized around 
0 and R = Ro^x, one should note that we have to anti¬ 
symmetrize two one-body states that are not orthogonal 
((0|R) 7 ^ 0), hence if 

|0,R)=AA(|0)®|R)-|R)®|0)) (B2) 

is the normalized and anti-symmetrized two-body state, 
(0, R|0, R) = A/'2 (2 (0|0) (R|R) - 2 (0|R) (R|0)) 

= 2A/'2 (l- |(0|R) 


(B3) 


the last term is not zero due to the non-orthogonality of 
coherent states, and it reads 


( 0 |R) = J d^r (r) 4 'r (r) = e ■» . 

Hence the normalization factor Af is 

- 1/2 


AA=(l-e-^)- 


(B4) 


(B5) 


such that the two-body state centered around 0 and R 
finally reads 




(ri,r 2 ) =A/'det 

\^o{r2) 4'R(r2)y 


27T 


0 4 0 2 ^0 


(r2-Ry 


(B 6 ) 


— 0 4 0 2^1^0 


(ri -R)^ 


The probabilty density of this state is given by 
p{r) = (i/)'*' (r)y>(r)) = 2 J \'i >2 (r,r 2 )|^ dr 2 


2tt 


(r-R)^ 


_r£ (r-R)^ ^ __ 

e 2 + 0 2 - 2 cos I I e 4 0 4 


r _x|-K^ 2 ^-Kg )2 r /y 2 Ro\ j 
X e 4 dx 2 / cos ( j e ^ dy 2 

(B7) 

This density can be written as the sum of two terms, 
pCr) = Po X) + 6p{r), with 


In order to compute the interaction energy of this state in 
the Hartree-Fock approximation we express this density 
in the reciprocal space. 


(BIO) 


Po ('^) = yy Po (!■) e dr 
= e“^ [l -I- e*'^"''^'’] 

/ Q^Rr\\ 

= 2M^ e- 5 — e"- cos f j 


and 


5 p(q) = / / 5p{v) dxdy 


= -2A/’^e g 2 ° Q \ cosh ^ ^y^o ^ ^ 

(Bll) 

This result can be generalized to a two-body state cen¬ 
tered around any pair of locations Ri and R 2 , 




(r„r2 ) = :^detf'^^^("i) ^Ir, (ri)\ 

ydlRi (l"2) 'I'R2(r2)y 


M 2 

X227, 


gjRiAri- 


(B12) 


_ ^^R2Ari- 


(*^1-^2) ip> (r 2 -Ri)^ 

4 4 


where 5R = ||Ri — R 2 II is the distance between the two 
centers and 


V 2 = (1 - e 2 j 


(B13) 


Finally 


P(r) = ^^ 


Ml 

2tt 






— 2 cos I -r A (5R • e- 


(B14) 


(r-R)^ _ 3 R^ 

2 0 8 


where R = (Ri -I- R 2 ) /2 et dR = Ri — R 2 . This density 
if shown on Fig. 18 for a two-body state localized around 
Ri = 0 and R 2 = 263 ,. As one might have expected, the 
main effect of the overlap correction is to suppress the 
electronic density between the two centers (i.e. at r = 
e^;), due to the antisymmetrization of the wave function. 

The corresponding Fourier transform is given by p = 
Po + (5p where 


M'^ 


e 2 -I- e 


(r-Hpe^c) 


(B 8 ) 


Po(q)=V|e-A ^ 

i^l,2 


(B15) 
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FIG. 16. Density without corrections 


2. Many-body generalization 
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The 7V-electrons Wigner crystal wave function is given 
by the Slater determinant of coherent states 

(ri,i’ 2 ,...,rAr) = A/^ det (^'r, (r^)) 

N 

= J\fN sgn(p) (rp(i)) , 

P£(7n 

where Ri are the N sites of the underlying lattice, and 
Afjsf is a normalization factor. The corresponding density 
is 


X [Ib] 


Pn (ri) = iV! 


I^'at (ri,r2,...,rjv)|^ dr2 ... dr^ 


FIG. 17. Density with corrections 



FIG. 18. (Color online) Corrections to the Wigner crystal 
density due to the overlap of neighbor states, localized around 
Ri = 0 and R 2 = 26^;. The top figure shows the uncorrected 
density for comparison. 


and 


= JVn X N\ X sgn(p) X sgn (p') 

P,p'6o’n 



(B17) 

At this step it is very common to suppose that the fol¬ 
lowing terms 


J (d) ^'r ,3 (!■*) dri oc 6a,p 


are negligible if a 7 ^ /3, i.e. the overlap of neighbor¬ 
ing localized states Rq and R /3 is neglected because 
their probability density decreases exponentially with 
distance. Within this approximation, which is justified 
in the low-density limit {ly ^ 1), the density (B17) is 
simply given by 


PN (ri) = En/■/ I^'r, dr 2 ... dr AT 

N 

= ^|vI/R. (r)^ 


2=1 


This density is the one of a direct product state 


(BI 8 ) 


5 p(q) = -Ee-^ JJ due-‘‘id"+R-) 

X cos ( E A 5R I e“ ^ 

V2 J 

= — 2 A/Je““’ *^e“^ e~^ cosh^^qA^R^. 

(B16) 

This correction is thus proportional to e“^*^ and gen¬ 
erally negligible compared to the main term poi when the 
two centers Ri et R 2 are separated by several magnetic 
lengths, because 


prod, (ri,..., Tat) = Tir^ (ri) T-r^ (r 2 )... T'r„ (rjv), 

which has a large overlap with the real fermionic 
state (B17) although it is not antisymmetric. 

For large densities, v ~ 1/2, the approximation to 
neglect the overlap of neighboring states is no longer 
valid. This is particularly true for the anisotropic crys¬ 
tals such as the rectangular and rhombic Wigner crys¬ 
tals, where neighboring sites may be very close to each 
other {SR ~ 21b)- While it would be tedious to take 
into account all the terms in the density (B17), we can 
at least take into account the overlap between nearest 
neighbors. Explicitely, under this approximation we sup¬ 
pose the terms 


<^P(q)/Po(q) oc e 


/ 




R„ 


(g) ^'Rp (d) dri 
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are non-zero if a = /3, as usual, but also if Re and R,g 
are nearest-neighbor sites of the crystal lattice. The in¬ 
tegrals involved in the calculation of the latter terms can 
be deduced from the two-body calculations presented in 
the previous section. 


The Fourier transform of the many-body state reads 


we obtain 


2e“V 

(5p(q) =-cosh 


1 o 

1 — 2 e 2 


(^)E 


iq R, 


1 o 

1 — 2 e 2 


e 2 e 2 cosh 


(^) 


X ^ ' 5 {Qx Qx) 'y ' ^ {Qy Qy) 

Qx Qy 

where Qx and Qy are defined by 


(B22) 


p(q) = 


1 o 

1 — 2 e 2 


e 2 


^(5(q- Qj) 


-2^ e-^e-^'i-^coshf ^qAJR 

(m^n) 


(B19) 


and 


27r 27r 

= V —n = —n 

6X v^i?o 


n 

Qy = = 


2'K^fa 

Ro 


m. 


(B23) 


(B24) 


where {m, n) refers to first neighboring sites, with relative 
coordinates (5R = Rm — Rn and center of mass R = 
(R„i -I- R„) /2, where Qj are the vectors of the reciprocal 
lattice. The normalization term [1 — 2exp(—(5i?^/2)]“^ 
was determined a posteriori with the constraint 


If we also take into account the site paired in the x 
direction, the correction to the projected density is 

5p (q) = - ^ 5 (q - Q) 

Q 


X 

-f 


4 - 2 e- 
2 


l- 2 e- 


5X2 


e ^2 cosh 


e cosh 



(B25) 


//x'> 


dr = p(q = 0) = iV, 


(B20) 


and the squared density 

l(p(q))l^ = l(p(q))l^ = e-«"/2^(I(q- Q) 

Q 

(B26) 

appearing in the energy thus becomes 


and we have used e“ = 

(27r)^ d (q — Qj )where ^„.c. is the unit¬ 
cell area of the Wigner crystal. 


We only consider a rectangular Wigner crystal that we 
have identified to be the crystal symmetry of lowest en¬ 
ergy. In a first step, we consider the pairs along the y 
direction, for which the lattice spacing 5R = 6YGy = 
R^Gy/y/a is decreased, with our choice a > 1 for the 
lattice anisotropy parameter. (Similarly, we have an en¬ 
hanced lattice spacing 5X = \/aRQ in the a;-direction.) 
Hence if we write the 


R^-fc = jSXGx + kSYGy {j, k} e (B21) 


l(p(q))l^ = 


1 - 2 


_ 

0 2+02 


)] '^^^(q-Q) 
Q 


- 9 V 2 


X 


1-4 



cosh 



-be 2 cosh 



-b4 

-b8 


^-SY- 


e 2 


^2fSY , 

cosh I —Qx ) -b e 

cosh (^^Qx ) cosh 


cosh^ 




-Qy 



(B27) 

The expression (B27) takes into account first-order cor¬ 
rections to the non-overlapping approximation (that is 
usually made when computing WCs energies). Glob¬ 
ally, these corrections energetically penalize strong 
anisotropies and/or large density, and we have used the 
expression (B27) for all our calculations. 
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